In an Ising model with spin-exchange dynamics damage always spreads 
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We investigate the spreading of damage in Ising models with Kawasaki spin-exchange dynamics 
which conserves the magnetization. We first modify a recent master equation approach to account for 
dynamic rules involving more than a single site. We then derive an effective-field theory for damage 
spreading in Ising models with Kawasaki spin-exchange dynamics and solve it for a two-dimensional 
model on a honeycomb lattice. In contrast to the cases of Glauber or heat-bath dynamics, we 
find that the damage always spreads and never heals. In the long-time limit the average Hamming 
distance approaches that of two uncorrelated systems. These results are verified by Monte-Carlo 
simulations. 
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I. INTRODUCTION 

Damage spreading (DS) investigates how a small per- 
turbation in a cooperative system changes during the 
time evolution (for a short review see, e.g., Ref. ||). 
In order to study DS two replicas of the system are con- 
sidered which evolve stochastically under the same noise 
realization (i.e. the same random numbers are used in 
a Monte-Carlo procedure). The difference in the micro- 
scopic configurations of the two replicas constitutes the 
" damage" . Depending on the Hamiltonian, the dynamic 
rules, and the external parameters the damage will either 
spread or heal with time (or remain in a finite spatial 
region). This behavior distinguishes chaotic or regular 
phases. 

Kinetic Ising models are among those systems for 
which DS has been studied most intensive. The ma- 
jority of the work has been devoted to single-spin-flip 
dynamic rules like Glauber, Metropolis or heat-bath dy- 
namics p|dil[ but also the Swendson-Wang cluster algo- 
rithm has been investigated jl2| . It has been found that 
the properties of DS (e.g. the question whether the dam- 
age spreads or heals for a particular model) depend sen- 
sitively on the dynamic rule chosen, i.e. DS is uniquely 
defined only if one specifies the Hamiltonian and the dy- 
namics. (Note that by considering all possible dynamic 
rules which are consistent with physics of a single replica 
an unambiguous definition of DS for a particular model 
can be obtained |ft3f.) 

The Glauber, Metropolis or heat-bath algorithms (as 
well as all other single-spin-flip algorithms) are exam- 
ples for a dynamics with non-conserved order parameter. 
There are, however, many physical systems that can be 
described by kinetic Ising models with order parameter 
conservation. A prominent example are, e.g., localized 
electrons where the Ising variables describe the electronic 
occupation numbers, and the dynamics consists of ther- 
mally assisted hops of an electron from one site to an- 
other. The simplest order parameter conserving dynam- 
ics in an Ising model is the spin-exchange dynamics of 
Kawasaki |L4|]. In this paper we want to investigate DS 
for this dynamics. To this end we first generalize the 



master equation approach ]10|Jll| to dynamic rules in- 
volving more than one site. We then derive an effective- 
field theory for DS in an Ising model with spin-exchange 
dynamics and solve it for a two-dimensional model on a 
honeycomb lattice. We find that in this model the dam- 
age always spreads. The stationary value of the damage 
is given by D* = (1 — m 2 )/2 (m is the magnetization) 
which corresponds to completely uncorrelated configura- 
tions. The results of the effective- field theory are con- 
firmed by Monte-Carlo simulations. 



II. MASTER EQUATION APPROACH 

We consider two identical Ising models with N sites 
described by the Hamiltonians and given by 
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where is an Ising variable with the values ±1, and 
n = 1,2 distinguishes the two replicas. is the ex- 
change interaction between the spins which we take to 
be J for nearest neighbor sites and zero otherwise. The 
dynamics (also called the Kawasaki dynamics |L4j) con- 
sists of exchanging spins on nearest-neighbor sites if the 
probability 



P = v(AE/2) 
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is larger than a random number £ G [0,1). Here AE 
is the energy change due to the exchange of the spins 
and T denotes the temperature. With this dynamics the 
total magnetization does not change with time, i.e. it is 
a conserved quantity. 

Within the master equation approach |lO|,[ll] the simul- 
taneous time evolution of the two replicas is described by 
the probability distribution 
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TABLE I. Damage creating and destructing processes for 
spin-exchange dynamics, all other processes do not change the 
damage. 
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where (•} denotes the average over the noise realizations. 
The variable Vi with the values H — h,H — , — h, or de- 
scribes the states of the spin pair (Sf,Sf). In the case 
of a spin-exchange dynamics the distribution P fulfills 
the master equation 
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P(vi, ...,u N ,t) = 

(ij) fJ-i.H 
(ij) W,« 
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where {ij} denotes all pairs of nearest neighbors and 
w(ui,Pj — > /ij, jij) is the probability for a transition of 
the states of the sites i and j from , Vj to fa, Uj. These 
transition probabilities can be obtained from (0). In ta- 
ble | we list all processes {vi, Vj — > fa, /ij) which lead to 
creation or destruction of damage, the probabilities for 
these processes will show up in the damage equation of 
motion later on. An important observation is that dam- 
aged sites can be created and destroyed only in pairs. All 
damage creating processes in table [I] can be transformed 
into each other by exchanging systems 1 and 2 and sites 
i and j . Their transition probabilities are therefore also 
related by symmetry. The same is true for all damage de- 
stroying processes. Thus, it is sufficient to calculate only 
two independent of the probabilities w(yi,Vj — > fa,fij), 
e.g., 
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where hi = J2j Jij^j 1S ^ ne l° ca l magnetic field of site i. 

As in the case of Glauber or heat-bath dynamics 
we derive an effective-field theory by assuming that 



fluctuations at different sites are statistically indepen- 
dent which amounts to approximating the distribution 
P(vi, ...,Vff,t) by a product of single-site distributions 
P Vi (t) . Order parameter conservation in the two systems 

imposes two conditions: P ++ (t) + P_| (t) = const and 

const. Inserting the decomposition 



P++(t) + P-+(t) 
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into the master equation (Q) gives a system of coupled 
equations of motion for the single-site distributions 
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where W(ui,Vj — * fa,fij) is the transition probability w 
averaged over the states v of all sites except for i and 
j. The total damage (Hamming distance) D can be ex- 
pressed in terms of the single-site distribution P v : 
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We note, that in contrast to Ising models with a non- 
conserved order parameter, the effective-field theory (Q) 
is not very useful in describing a single system since the 
only remaining dynamic variable for a single system, viz 
to, does not change during the time evolution. The dam- 
age is, however, not conserved and (Q) constitutes a useful 
mean-field theory for its time evolution. 

From the single-site master equation ([?]) and table | we 
derive an equation of motion of the damage. Using some 
symmetry relations |l5| between the transition probabil- 
ities W, it reads 



— D = -2D 2 W(+-, 
+2[(1 - D) 2 - to 2 ] W{++, 



++,—) (9) 



So far the considerations have been valid for all di- 
mensions and lattice types. To proceed we now study 
a particular system, viz. a two-dimensional Ising model 
on a honeycomb lattice. The same model was studied 
in Refs. jl0|,[llj for single-spin-flip dynamic rules. The 
calculation of the average transition probabilities W can 
be carried out analogously to the case of single-spin-flip 
dynamics. It is straightforward but tedious, the details 
will be published elsewhere fljsf . Here we only give the 
results: 
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Here t„ = tanh(n J/T). If we insert (10b) into the dam- 
age equation of motion (^|) we observe that the death 
term [first line of eq. (||)] is of order D 2 for small D. This 
is a major difference to the case of Glauber and heat bath 
dynamics [fjoj[ll| (where the death term is of order D) 
and reflects the fact that for spin-exchange dynamics the 
damage can only be destroyed pairwise. In contrast, the 
birth term is of order D (as it is for Glauber oder heat 
bath dynamics) because already a single damaged site 
can produce further damage in its neighborhood. Con- 
sequently, for small enough D the birth term will always 
be larger than the death term and the damage will never 
heal completely. 

We now discuss the stationary solutions of the damage 
equation of motion (||) and their stability. We restrict 
ourselves to the case that the system is in equilibrium 
when the damage is introduced. Thus, to can be taken 
to be the equilibrium value of the magnetization which 
is zero for T > T c « 2.11J and 
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for T < T c in our effective- field theory frL0]j . Obviously, 
D = D{ = is always a fixed point (FP) of (§). To 
investigate its stability we expand (|]) to linear order in 
D. The resulting linearized equation of motion is given 
by dD/dt = X X D with 

Ai = (1 - to 4 ) tanh ^ + i(l - to 2 ) 2 tanh (12) 

The Lyapunov exponent Ai is always positive, thus the 
FP D\ is always unstable. In Fig. [j] we show the temper- 
ature dependence of Ai. Since the Hamiltonian and the 
dynamic rule are invariant under a global flip of all spins, 
the existence of the FP D\ = implies the existence of 
the FP D\ — 1 with the same stability properties. 

In the paramagnetic phase the only other stationary 
solution of @ in the physical interval < D < 1 is 
Z?3 = 1/2. By expanding (||) around D = 1/2 we obtain 
the corresponding Lyapunov exponent, 
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Since A3 is always negative the FP D3 = 1/2 is stable in 
the entire paramagnetic phase. The temperature depen- 
dence of A3 is also shown in Fig. |]. 

In the ferromagnetic phase there are two more sta- 
tionary solutions of (0) in addition to D* = and 
£>2 = 1- If the two systems have the same value of 
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FIG. 1. Results of the effective-field theory for damage 
spreading in the kinetic Ising model with spin exchange dy- 
namics. 

the magnetization, to/ 1 ' = to/ 2 ) = to we obtain the 
FP D3 = (1 — to 2 )/2. If the two systems have oppo- 
site magnetization to/ 1 ) = —m^ = to we obtain the 
FP D\ = (1 + to 2 )/2. Since D* and D\ are related by a 
global flip of all spins in one of the systems, they have the 
same stability properties. The corresponding Lyapunov 
exponents are given by 
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They are always smaller than zero, thus and D\ are 
stable in the entire ferromagnetic phase. The tempera- 
ture dependence of D3 and A3 is shown in Fig. |l} 

We now show that at the stable FPs the configurations 
of the two systems are completely uncorrelated. From the 
definition (M) of the damage we obtain 
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For uncorrelated configurations of and we have 

(S^S^) = (3^ l) )(Sl 2) ) = ±to 2 if the two systems have 
equal or opposite magnetization, respectively. Thus, for 
uncorrelated configurations we obtain D = (1 =F to 2 )/2. 
These are exactly the values of D| and D\. 



III. MONTE-CARLO SIMULATIONS 

We have verified the main predictions of the mean- 
field theory by Monte-Carlo (MC) simulations of a three- 
dimensional Ising model with Kawasaki spin-exchange 
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FIG. 2. Time evolution of the damage in an Ising model 
with spin-exchange dynamics. The data points represent av- 
erages over 100 runs of a system of 27 3 sites. 



dynamics according to (|^). The simulations are carried 
out on cubic lattices with up to N = 101 3 sites with pe- 
riodic boundary conditions. By comparing different sys- 
tem sizes we verify that any finite size corrections to the 
results are smaller than the statistical error of the sim- 
ulation. This is easily possible since we are away from 
a spreading transition and thus the damage correlation 
length is finite. 

In this study we are not interested in phase separa- 
tion processes. We thus prepare the system with the 
correct equilibrium magnetization value for each temper- 
ature simulated. 2000 MC sweeps are carried out to equi- 
librate the system. Then the initial damage Dq is created 
by exchanging randomly chosen pairs of nearest neighbor 
spins in one of the systems. We use values of Do between 
5 * 10~ 4 and 5 * 10~ 2 . After that both systems evolve in 
parallel using the same random numbers. Examples of 
the time evolution of the damage are shown in Fig. |^. 
Within the first 5 to 10 MC sweeps the damage increases 
approximately exponentially with time. A fit of the data 
to an exponential law gives an estimate for the Lyapunov 
exponent Ai. We note that the damage time evolution 
shows a systematic deviation from an exponential law 
which manifests in a slight downward curvature in Fig. 
|2|. This deviation stems from the fact that we are simulat- 
ing a lattice system. Since with Kawasaki dynamics the 
damage can spread at most two lattice constants per time 
step the increase of the damage with time is bounded by 
a power law, D(t) < D max ~ (2t) 3 Do. Therefore, a pure 
exponential spreading can only be observed as long as 
the probability for any site (or pair of sites) to become 
damaged during a particular time step is small compared 
to one. (In this case the above bound set by the lattice 
does not play a role.) For our system this condition is, 
however, only fulfilled for small temperatures. 

In order to determine the long-time limit of the aver- 
age damage we average its values over 5000 MC sweeps 




FIG. 3. Asymptotic average damage DJ and Lyapunov ex- 
ponent Ai for the kinetic Ising model with spin exchange dy- 
namics. The Lyapunov exponents have been obtained from 
100 runs of a 27 3 system. The FP values _D| has been calcu- 
lated from 10 runs of a 101 3 system. Their statistical errors 
are smaller than the symbol size. 



after a plateau has been reached. The results of our sim- 
ulations are summarized in Fig. 3| We indeed find that 
the FP D\ is unstable, and the damage always spreads. 
The Lyapunov exponent Ai of the FP D\ is positive for 
all temperatures investigated. The asymptotic average 
damage takes exactly the value of two uncorrelated con- 



m 2 )/2. (We did not observe 
(1 + m 2 )/2 since we always 
started with the two systems having the same magneti- 
zation.) 



figurations, viz. D\ = (1 
the other stable FP D\ - 



IV. CONCLUSIONS 

To summarize, we have used an effective-field theory 
and Monte-Carlo simulations to show that the time evo- 
lution of a kinetic Ising model with Kawasaki spin ex- 
change dynamics is chaotic for all temperatures in the 
sense that the FP D\ = is unstable. Moreover, we have 
shown that two systems whose initial configurations dif- 
fer only at a few sites become completely uncorrelated in 
the long-time limit. This corresponds to an asymptotic 
average damage of D = (1 — m 2 )/2. 

In this last part of the paper we want to discuss how 
general these results are. Since the properties of DS are 
known to depend on how the random numbers are used 
in the update process |l3(] for single-spin-flip dynamics, 
an analogous comparison for spin-exchange dynamics is 
desirable. However, the main properties of our solution 
will be robust against such changes in the update rules. 
In particular, the fact that the damage death rate (see 
eq. ||) is of order D 2 is a result of the spin-exchange 
mechanism alone. It is therefore independent of how the 
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FIG. 4. Comparison of the asymptotic average damage D't, 
for the kinetic Ising model with Kawasaki and modified spin 
exchange dynamics. 



random numbers are used in the update rule. This sug- 
gests that the main finding of this paper, viz. that the 
spin-exchange dynamics is chaotic for all temperatures is 
valid not only for the Kawasaki update rule (0) but in 
general. As a first step of a future systematic investiga- 
tion of different update rules we have studied a modified 
version of the Kawasaki dynamics. The modification con- 
sists of using the random number £ if the configuration 
of the spin pair selected for the exchange is (H — ) but 
using 1 — £ instead if the configuration is ( — h). This 
modified update rule can be seen as the spin-exchange 
analog of the heat-bath dynamics (in the same sense as 
the Kawasaki dynamics can be seen as the analog of the 
Glauber dynamics) . In Fig. |] we compare the asymp- 
totic average damage of the Kawasaki and the modified 
spin-exchange dynamics. The modified dynamics gives 
lower damage values than the Kawasaki dynamics for all 
temperatures. Nonetheless, the fixed point D\ = is 
unstable for all finite temperatures, and the asymptotic 
damage is finite. This is in agreement with the above sug- 
gestion that a spin-exchange dynamics is always chaotic 
irrespective of the particular update rule. 

Let us finally discuss the relation of the DS process dis- 
cussed here with other non-equilibrium processes. As al- 
ready mentioned, a key feature of DS with spin-exchange 
dynamics is that damaged sites can heal only in pairs 
while they can diffuse alone and also create further dam- 
age. This is different from the contact process and other 
processes in the directed percolation universality class 
where a single active site can die locally with finite prob- 
ability. There is, however, a simple reaction-diffusion 
process which should show qualitatively the same behav- 
ior as DS with spin-exchange dynamics. Since for small 
damage D the birth rate is proportional to D (and since 
damage is created in pairs) while the death rate is pro- 
portional to D 2 (see eqs. (P,[l0|), such a reaction-diffusion 
process could be defined by the reactions 



and additional diffusion of the substance A. For small 
concentrations of A it should have the same qualitative 
behavior as DS with spin-exchange dynamics for small 
damage. 
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